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ABSTRACT 

The quasar remnant cores of nearby giant elliptical galaxies NGC 4486 (M87), 
NGC 1399, NGC 4649 and NGC 4472 are the sites of supermassive (> 10 9 
Mq) black holes. These objects are investigated as to the viability of the 
conjecture that they could harbor compact dynamos capable of generating the 
highest energy cosmic rays. For an accretion process involving an equipartition 
magnetic field near the event horizons of the underlying putative spun-up 
black holes, the energy achievable in accelerating protons could well be > 
10 20 eV for all these when only considering the drag induced by curvature 
radiation. Estimates of the SED (spectral energy distribution) of ambient 
core photons lead to the conclusion that the energy losses arising from photo- 
pion production in proton collisions with these target photons are relatively 
small for all but M87. For M87, the ambient photon field is likely to be a 
limiting factor. Accretion rates of ~ 1 Mq yr _1 , comparable to the Bondi 
rates and to the stellar mass loss rates, are associated with (> 10 20 eV) cosmic 
ray generation in the other (electromagnetically dark) galactic core sites. If 
these sites are found to be sources of such cosmic rays, it would suggest the 
presence of a global inflow of interstellar gas all the way into the center of the 
host galaxy. 

Key words: acceleration of particles - black hole physics - cosmic rays - 
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1 INTRODUCTION 

The massive dark objects (MDOs) at the centers of present-epoch giant elliptical galaxies 
appear to be the supermassive black hole remnants of earlier quasar activity (McLure et al. 
1999; McLeod, Rieke & Storrie-Lombardi 1999; Salucci et al. 1999; Richstone et al. 1998; 
Boldt &; Leiter 1995; Chokshi & Turner 1992). It has been proposed that the highest energy 
cosmic ray particles, those above the 'GZK' limit (Greisen 1966; Zatsepen & Kuz'min 1966), 
could be accelerated by the effective electromotive force (emf) generated near the event 
horizon of spinning supermassive (> 10 9 Mq) black holes at the apparently dormant quasar 
remnant cores of nearby massive elliptical galaxies QBoldt fc Ghosh 1999| ). While an emf up to 
10 21 volts or somewhat greater would appear to be possible, the drag arising from curvature 



radiation ( |Levinson 2000| ) limits the energy to be attained by a charged particle accelerated 
by this field; for the four cases considered here (i.e., NGC 4486, NGC 1399, NGC 4649 & 
NGC 4472) the suppression factor for a proton is estimated to be ~ 6 — 9. Acceleration 
to energies > 10 20 eV would then still be feasible. However, during the final phase of the 
acceleration process energy losses due to inelastic collisions with ambient photons, those 
dominated by photo-pion production, must be considered. Therefore, we address in some 
detail the major concern that recent radio data on low-luminosity AGNs could imply a 
photon number density at the nuclear core that is so high as to preclude ultra-relativistic 
processes of particle acceleration and escape, particularly as they pertain to the immediate 
vicinity of the putative supermassive black hole flBlandford 200C| ). 

Since 7 = 10 11 for a 10 20 eV proton, ambient photons at hv > 15 x 10~ 4 eV (i.e., 
v > 360 GHz; A < 0.83 mm) would appear in the nucleon's proper frame as 7-rays > 150 
MeV, above the critical threshold for energy losses via inelastic collisions involving pion 
production (cf., Stecker 1968; Hill & Schramm 1985; Stecker & Salamon 1999). Using recent 
estimates of the SED (spectral energy distribution) at the core of M87 (NGC 4486), the 
central (and dominant) giant elliptical galaxy in the Virgo Cluster (Reynolds et al. 1996; 
Ho 1999; Di Matteo et al. 2000), we find that the corresponding radiation length of such 
energetic protons might well be significantly smaller than the Schwarzschild radius, thereby 
precluding the acceleration to such a high energy. Hence, if M87 is indeed found to be a 
source of the highest energy cosmic rays, as commonly conjectured (Biermann & Strittmatter 
1987; Biermann 1999; Farrar & Piran 2000), the underlying accelerator would then have 
to be other than the compact dynamo considered here. However, for other nearby giant 
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elliptical galaxies harboring apparently dormant supermassive (> 10 9 Mq) black holes, in 
particular for the remaining three such objects addressed by Di Matteo et al. (1999, 2000), 
the radiation lengths are estimated to be larger than their Schwarzschild radii. These include 
NGC 1399, the central giant elliptical in the Fornax Cluster and two other giant elliptical 
galaxies in Virgo (i.e., NGC 4649 and NGC 4472). The dynamo characteristics expected for 
the associated compact cores (i.e., 5-field, achievable energy) are derived from independent 
estimates of the black hole mass, Bondi accretion rate and SED environment. It has been 
pointed out that TeV 7-ray curvature radiation is a necessary consequence of cosmic ray 
generation by such black hole dynamos ( [Levinson 2000| ); these giant elliptical galaxies would 
be worthwhile targets for the observation of such electromagnetic radiative signature. 



2 DYNAMO CHARACTERISTICS 

The dynamo's emf (V) is generated by the black hole induced rotation of externally supplied 
magnetic field lines threading the horizon ( [Blandford fc Znajek 1977|) . If B is the ordered 
poloidal field near the hole, V ~ aB, where a is the hole's specific angular momentum; for a 
hole mass M, a < M (e.g., a = M for an extreme Kerr hole). In astrophysical units flZnajek 
19781) : 

V = 9 x 10 20 (a/M)B 4 M 9 volts, (1) 

where B 4 = £/(10 4 G) and M 9 = M/(1O 9 M ). 

The energy density of the magnetic field near the event horizon is expected to be in 
equipartition with the rest mass energy density of accreting matter ( [Krolik 1999Q . In terms 
of an Advection Dominated Accretion Flow (ADAF) model (cf., Di Matteo et al. 1999) this 
is to be identified with the regime where the gas pressure is half the total (i.e, (3 = 1/2). 
Under this assumption, 

B A = l^M^M 1 / 2 , (2) 

where M is the accretion rate dM/dt in Mq yr _1 . 

The magnetic field (B) to be associated with the nuclei of those elliptical galaxies con- 
sidered here is obtained via equation 2 using the Bondi accretion rates Msondi estimated by 
Di Matteo et al. ( |2000|) and the black hole masses determined by Magorrian et al. ( |1998| ). 
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Noting that the hole's radius is fixed by its mass, the emf given by equation 1 may be 
expressed directly in terms of the accretion rate. The maximum emf (V), that corresponding 
to a/M close to unity, is then given by: 

V = 1.2x 10 21 M 1/2 volts. (3) 

Except where noted, we take M = M Bon di- 

The energy (E) to be attained in this electric field is limited by the drag arising from 
curvature radiation induced by the magnetic field ( |Levinson 2000| ). From equation 5 in 
Levinson (op. cit.), we obtain that, for a proton (charge e + ), the suppression ratio is given 
by: 

E/[e(V)} « [(50M 9 )- 1/2 B 4 " 3/ V /2 > (4) 

where r is the magnetic field curvature in units of the Schwarzschild radius. For r ~ 1 and 
M « (0.1 - 10) M yr" 1 , we note (from equations 2-4) that E = (1.0 - 1.8) x 10 20 M 9 1/4 
eV. Although the energy possible for a heavier nucleus could be substantially greater, such 
a particle (e.g., Fe) launched at > 10 21 eV would be disrupted into its constituent nucleons 
after traveling only 20 Mpc (e.g, Cronin 1997). By contrast, a proton starting with 10 21 eV 
would maintain about one-third of its initial energy after traversing the same distance. 

During the final phase of the acceleration process energy loss due to photo-pion produc- 
tion in collisions with ambient photons becomes a relatively important effect . The associated 
radiation length (A) relative to R (the radius of core emission) is given by 

A/R = C7rR/((Ka}Q) (5) 
= (278/(AV))(i^s)M 9 (Q/10 53 a" 1 )" 1 , 

where a is the proton photo-pion production cross-section, o' its value in microbarns (10 -30 
cm 2 ), Rs is the Schwarzschild radius, K = (E(loss)) /^(initial) is the inelasticity in a single 
collision ( Stecker 1968 ) and Q is the core emission rate (photons s -1 ) for electromagnetic 
radiation at v > 360 GHz, given by 

Q = h- 1 J v- x L v dv (6) 

where h is the Planck constant and L u = 4nD 2 E u for a source of spectral density F v at 
distance D. We note ( |Uaso et al. 1998 ) that, for the regime of interest here, 



(Ka>) 



(Ka'(dQ/du)du 



/Q < 120 microbarns. (7) 
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Table 1. Black Hole Galactic Nuclei: Candidate Cosmic Ray Sources 



Host Galaxy 
Host Cluster 




NGC 4486 
Virgo 


NGC 1399 
Fornax 


NGC 4649 
Virgo 


NGC 4472 
Virgo 


Distance (D) 


(Mpc) 


18 


29 


18 


18 


Black Hole Mass (AT) 


(10 9 M ) 


3.6 


5.2 


3.9 


2.6 


Meondi 


(M yr" 1 ) 


1.5 


3 


1.4 


0.7 


Stellar Bulge Mass (M Bulge ) 


(10 12 M ) 


0.82 


0.32 


0.54 


0.84 


-kobs/^Bondi 




lO" 5 


2 x 10- 6 


3 x 10- 5 


lO" 5 


Magnetic Field (B) 


(10 4 G) 


0.45 


0.44 


0.40 


0.43 


emf (V) 


(10 20 volts) 


15 


21 


14 


10 


r -V* E 


(10 20 eV) 


_ 


2.4 


2.0 


1.7 


Obs. Radio Freq. (f b s ) 


(GHz) 


100 


43 


43 


43 


VLt v (t'obs ) 


f"10 38 ertr s — 1 ) 


32 


6.7 


2.0 


3.3 


<i/L„(360 GHz)) 


(10 38 erg s" 1 ) 


47 


13 


3.8 


6.2 


vL u {7 x 10 8 GHz) 


(10 38 erg s" 1 ) 




< 13 






(Q(A < 0.83mm)) 


(10 53 photon s" 1 ) 


~ 28 


< 7.7 


< 2.3 


< 3.7 


(R S /R)A/R 




~ 0.3 


> 2 


> 4 


> 2 



The black hole and galactic bulge masses are from Magorrian c t al. ( 199? 1. The distances, Bondi accretion 



rates and ^obs/^Bondi ratios are taken from Di Matteo et al. (200C). Radio observations are at 0.3 cm for 
NGC 4486/M87 (Reynolds et al. 1996; Baath et al. 1992) and 0.7 cm for the three other sources (Di Matteo 
et al. 1999, 2000). X-ray observations with the Chandra Observatory at hu = 3 keV (Locwcnstcin 200C) 
correspond to v = 7 X 10 8 GHz. 

3 CONSTRAINTS 



The key characteristics of the four dynamo candidates investigated here are summarized in 
Table 1. The tabulated values for vL v at v = 360 GHz were obtained from F v measurements 
at lower frequencies under the assumption that F v oc v~ ' 7 , consistent with the radio spectra 
observed in the vicinity of 43 GHz (Di Matteo et al. 1999, 2000; Reynolds et al. 1996). The 
core emission rate (Q photons s" 1 ) at A < 0.83 mm {y > 360 GHz) from the accretion 
flow associated with each of the four supermassive black holes is estimated via equation 6 
under the assumption that F u oc z/~ ar at higher frequencies as well. For M87 this is based 
on a high resolution (10 -4 arcsec) VLBI measurement at 100 GHz (Reynolds et al. 1996; 
Baath et al. 1992). For the other three sources we rely on VLA measurements at 43 GHz (Di 
Matteo et al. 1999, 2000). Because these latter measurements allow more of a contribution 
from the underlying galaxy and weak jets, and two (NGC 4472 & NGC 4649) do show some 
extended emission, they are taken to be upper limits to the flux from the compact cores of 
these galaxies. Comparisons of these extrapolated spectra with the observed limits shown in 



© 2000 RAS, MNRAS 000, L00-L00 



L6 E. Boldt and M. Loewenstein 

Figure 2 of Di Matteo et al. ( |2000| ) are consistent with these estimates being upper limits to 
L v at v > 360 GHz, although less clear for M87. The presence of inner jets associated with 
any of these potential dynamos would imply vacuum breakdown if created by the mechanism 
suggested by Blandford & Znajek (1977). In this situation, the feasibility of our compact 
cosmic ray generator demands that such a jet ejection process be episodic. 

The magnetic field (B) is obtained from the Bondi accretion rate and black hole mass via 
equation 2. The emf is calculated from the Bondi accretion rate via equation 3. Considering 
the drag arising from curvature radiation flLevinson 2000|) the energy (E) to be attained by 
a proton accelerated by this emf is estimated via equation 4. During the final phase of the 
acceleration process energy losses due to collisions with ambient photons are dominated by 
photo-pion production. The associated radiation length (A) relative to R listed in Table 1 is 
obtained from equation 5, where R is the radius of core emission and R$ is the Schwarzschild 
radius. 

Lack of microwave data for these sources prevents any direct confirmation that our power- 
law extrapolation from the radio (centimeter) band does in fact provide an upper limit to 
their emission in the entire relevant submillimeter band at A < 0.83 mm. As a result, we 
depend on the ADAF models prescribed by Di Matteo et al. ( |2000| ) to provide us with 
the template needed for evaluating our power-law extrapolation. For M87 our extrapolation 
equals or exceeds L u expected for all their models except the 'no-wind' ADAF; see Figure 2a 
in their paper. For the three other sources our power-law definitely exceeds what is expected 
for all models except the 'no-wind' ADAF; see Figures 2b, 2d & 2e. Radio data exhibited 
by Di Matteo et al. (|2000l) at longer wavelengths, at least for NGC 4649 and NGC 4472, 
appear to rule out the no-wind ADAF model in favor of models with winds that have less 
emission at frequencies > 360 GHz than those based on the simple extrapolations presented 
in Table 1. Moreover, analysis of Chandra X-ray Observatory data for NGC 1399 (see below, 
and Loewenstein et al. 2000) places an upper limit on vL v at 3 keV (~ 7 x 10 17 Hz) of 
1.3 x 10 39 erg s _1 (where, for consistency, we adopt the distance to NGC 1399 of 29 Mpc 
used by Di Matteo et al. 2000) with comparable limits throughout the 0.3-3.0 keV X-ray 
energy range. A preliminary estimate of the detected UV flux from a nuclear point source 
yields a comparable value of vL v at v ~ 2 x 10 15 Hz ( P'Connell 2000| ). Simply extrapolating 
the highest energy radio point (at 43 GHz) through the X-ray upper limit reduces the value 
of the photon emission rate Q in Table 1 by about a factor of 3 for NGC 1399, thereby 
increasing the associated lower limit to A/R threefold. 
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4 IMPLICATIONS 

The vicinities of spun-up supermassive black holes in dormant elliptical galaxy nuclei can 
provide an acceleration mechanism sufficient to explain observations of the highest energy 
cosmic rays, and meet the further criterion of being numerous in the nearby (< 50 Mpc) 
universe defined by the GZK cutoff for the highest energy particles. Since few, if any, other 
astrophysical sites in the present-epoch are as feasible, the comparison of the following 
predictions with cosmic ray data of improved spatial and spectral statistics could deter- 
mine whether models invoking new particles or topological defects need be considered more 
carefully Conversely, confirmation of the present model has far-reaching astrophysical im- 
plications; e.g, it would indicate continued spin-up and growth of supermassive black holes 
via accretion. 

The angular distribution of ultra-high energy (UHE) cosmic ray sources provides the 
cleanest test of a supermassive black hole origin, as the most massive black holes lie in the 
nuclei of elliptical galaxies. Scattering by intracluster fields (microgauss B) may destroy 
this association for galaxies in clusters; however, TeV 7-rays created via curvature radiation 
( [Levinson 2000| ) would preserve the correlation. Also, one would not expect otherwise active 



galaxies, with their higher central photon densities, to produce these cosmic rays. Finally, 
under the equipartition assumption the emf should be proportional to the square root of 
the rate of accretion into the nucleus, which should be well approximated by the Bondi rate 
(Brighenti & Mathews 1999; Quataert & Narayan, 2000), i.e. V oc p g LT sas M, where p gas 
and T gas at the accretion radius are measurable using the Chandra X-ray Observatory. 

The estimates of the radiation lengths in Table 1 rest on the extrapolation of the radio 
spectrum through 360 GHz assuming a z/ -0 7 energy spectrum. A significantly flatter (or 
inverted) slope would imply a photon number density that is prohibitively high for particle 
acceleration to > 10 20 eV; and, indeed this seems to be the case for many low-luminosity 
AGN ([Ho 1999| ). However, with the exception of M87, the galaxies considered here may 



not be members of the same population, as they lack detectable optical emission lines ( Ho, 
Filipenko fc Sargent 1997| ). 



The Chandra limit on the 3 keV nuclear emission for NGC 1399 is derived using the 
extracted spectrum in a 3.5 arcsec diameter aperture. Two-component models including 
thermal (hot ISM) and power-law components are fit to the data and 90% confidence upper 
limits to the power-law flux derived. This upper limit is robust as it accounts for all the 
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observed emission above 2 keV, and is conservative since image analysis strongly suggests 
an even lower upper limit to the flux originating in any nuclear point source ( [Loewenstein 
20001 ). The hard component attributed by Di Matteo et al. ( |2000| ) to an ADAF associated 
with the nuclear supermassive black hole has now been resolved with Chandra into an 
extended distribution of discrete sources ( |Angelini et al. 2000| ). The VLA radio observations 
of NGC 4649, NGC 4472, and NGC 1399 correspond to a similar angular resolution as that 
of Chandra and should be regarded as upper limits as well since there may be contributions 
from inner-jet emission [the spectra in Ho (|1999|) are derived from higher resolution VLBI 
observations]. Thus these galaxies are quiescent even relative to LINERS and other low 
luminosity AGN, and it is plausible that cosmic ray energy losses by ambient photons are 
unimportant for these three galaxies. We note that M87, a likely exception, displays more 
powerful activity at all energies. 

If the Chandra results for NGC 1399 are generalizable, the nuclear X-ray emission from 
elliptical galaxies is well below that predicted by models of ADAFs accreting at the Bondi 
rate. However, for the specific cases considered in this paper the Bondi accretion rate that 
provides sufficient compression of the magnetic field to generate an emf corresponding to 
the most energetic cosmic rays is of the same order as the integrated stellar mass loss rate. 
M B ondi/M B uige = 1-8, 9.4, 2.6 and 0.83 xl0~ 12 yr" 1 for M87, NGC 1399, NGC 4649 and 
NGC 4472, respectively (Table 1). These are comparable to the specific mass return rate 
calculated for the stellar population in elliptical galaxies (e.g., Mathews 1989), suggesting 
the presence of a global inflow of interstellar gas that persists from large galactic radii all 
the way into the very center of the galaxy. Thus while (hot) protogalactic gas may be the 
source for the initial rapid growth and quasar-epoch fueling of massive black holes ( Nulscn 
fc Fabian 2000| ), the stellar population of the galaxy at large may be responsible for 'feeding 
the monster' ( |Gunn 1979|) during the present era: an era where the quieter power output may 
very well be characterized by high energy particles and their associated gamma radiation. 



5 OUTLOOK 

Comprehensive investigations of nearby giant elliptical (& SO) galaxies indicate that most 
harbor a massive dark object (MDO) at their centers ( |Magorrian et al. 1998j ). Studies of 
AGN evolution (Chokshi & Turner 1992; Fabian & Iwasawa 1999; Salucci et al. 1999) and the 
X-ray background of accretion-powered radiation flBoldt fc Leiter 1995|) conclude that the 
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largest of these MDOs are associated with supermassive black holes which are the apparently 
dormant remnants of previously active quasars. It has been proposed that these accretion- 
fed spinning dark objects are latent dynamos, sufficient for producing the highest energy 
cosmic rays ( |Boldt fc Ghosh 1999| ). The accretion rate for such a dynamo is no greater than 
the mass loss rate estimated for giant elliptical galaxies ( |Mathews 1989[ ). The scatter and 
uncertainty in the ratios of black hole to bulge mass and of bulge mass to optical luminosity, 
and the inhomogeneity of the galaxy distribution within 50 Mpc, make it virtually impossible 
to precisely quantify the expected underlying number of UHE cosmic ray sources. However, 
we note that the number of galaxies within 50 Mpc that have bulges sufficiently luminous 
(-^Buige > 10 10 L b q) to be potential sites of > 3 x 10 s Mq black holes is on the order of 10 3 
(Magorrian 1999, Marinoni et al. 1999). Assuming an IGM B field of coherence length ~ 1 
Mpc and magnitude ~ 1 nanogauss, a UHE proton > 10 20 eV originating within 50 Mpc 
would have an arrival direction aligned sufficiently well with its origin for avoiding confusion 
among candidate sources QMedina-Tanco, de Gouveia Dal Pino, fc Horvatth 2000| ) , although 
it has been suggested that the relevant B field could be larger than a nanogauss flFarrar 
fc Piran 2000| ). For a UHE proton < 10 21 eV originating from a source deep within a rich 
cluster of galaxies the Larmor radius (for the associated microgauss B field) is much less than 
the cluster radius QSaikia fc Salter 1988|) ; hence the emerging protons would appear to be 
coming from an extended source the size of the cluster. Fortunately, most suitably massive 
giant elliptical galaxies reside outside of clusters (Magorrian et al. 1998; Burstein 1999). If the 
'local' extragalactic magnetic field is much less than a microgauss, then an angular resolution 
of about a degree should be good enough for the OWL/Airwatch air-shower observatory 
( [Strcitmatter 1998 ) to establish a correlation between candidate ellipticals and UHE cosmic 
ray events, at least two-thirds of the time (the remainder would be cluster associated). Those 
in clusters would then be located via their TeV 7-ray curvature radiation. In particular, of the 
six such giant ellipticals addressed by Di Matteo et al. (|2000|) we find that, after considering 
principal constraints, at least half could well be harboring supermassive black hole dynamos 
capable of producing cosmic rays more energetic than 10 20 eV and be identifiable via their 
TeV 7-radiation, e.g., with the Whipple Observatory at an angular resolution sufficient for 
isolating specific galaxy emission ( |Weekes et al. 199~6| ). 
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